Spin ordering quantum transitions of superconductors in a magnetic field 
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We argue that recent neutron scattering measurements by Lake et al. (Science 291, 1759 (2001)) 
of the spin excitation spectrum of La2-5Sr5Cu04 in a magnetic field can be understood in terms of 
proximity to a phase with co-existing superconductivity and spin density wave order. We present a 
general theory for such quantum transitions, and argue that their low energy spin fluctuations are 
controlled by a singular correction from the superflow kinetic energy, acting in the region outside 
the vortex cores. We propose numerous experimental tests of our theory. 
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Recent neutron scattering experiments of Lake et al. |^] 
have opened a new window on the spin excitation spec- 
trum of the high temperature superconductors. They 
have observed that a moderate magnetic field induces 
an anomalously large increase in the spectral density of 
low energy spin fluctuations in La2-5Sr5Cu04 at opti- 
mal doping {6 = 0.163) and low temperatures (T). Ex- 
periments on the underdoped La2-5Sr5Cu04 also show a 
large increase in the intensity of elastic neutron scatter- 
ing in the applied magnetic field There have been 
a number of studies of enhanced antiferromagnetic cor- 
relations in the cores of vortices in the superconducting 
order |^-^ , and one interpretation of the measurements 
is that the extra scattering arises from quasi-static mo- 
ments in the cores of the field- induced vortices. The 
measurements then lead to the estimate [Q that the each 
Cu site in the vortex core has a moment of order 0.6/is 
in the anomalous low energy sector. Such a large mo- 
ment is characteristic of the insulating antiferromagnet 
{S = 0), and would require a corresponding charge dis- 
proportionation with a large Coulomb energy cost; this 
is a difficulty with such an interpretation. 

We argue here that the experiments can be under- 
stood by assuming that the superconductor (SC) is in 
the vicinity of a bulk quantum phase transition to a 
state with microscopic co- existence of SC and spin den- 
sity wave (SDW) orders; the latter state has been consid- 
ered in a number of studies [§-0 , and has been observed 
in excess-oxygen-doped La2Cu04+y jl^. The magnetic 
field, H, drives the SC phase closer to the SC-hSDW 
phase — see Fig ^ Initially, H induces well separated vor- 
tices in the SC, and there could be small additional mag- 
netic scattering from relatively high energy spin S = 1 
excitons centered around the vortex cores. However, with 
increasing H (but with H/H^2 still small the up- 

per critical field at which the mixed vortex state dis- 
appears for a particular set of parameters-see Fig |^)) 
the energy of any such exciton decreases rapidly, and it 
becomes part of a delocalized collective spin fluctuation 
which is a precursor of the transition to the SC-I-SDW 
phase; the dominant magnetic scattering then arises from 



the region outside the vortex cores. Our primary results 
for extreme type-II superconductors (Ginzburg-Landau 
parameter k oo) are: (i) for small HlH%^, the SC 
to SC-I-SDW phase boundary (MA in Fig y the caption 
defines r, Tc) is at iJ ^ (r — rc)/ln(l/(r — r^) — so this 
boundary is anomalously flat at small H , and this allows 
the system to move close to it for not too large H; (zi) 
the energy, eoo(-ff), of the low energy peak in the neu- 
tron scattering cross section should decrease as cqo (i? ) = 
eoo(O) - Ci{H/H°2)ln{H^2/H) along the vertical arrow 
in Fig |l|, (Hi) when starting from the SC-I-SDW phase 
at H = 0, one finds an increase of the elastic neutron 
scattering given by /(ff) = I{0)+C'2{H/ H^2)HHc2/ H), 
where Ci and C2 are some constants. All these functional 
forms are expected to be exact at small H; naively, one 
might have expected an analytic series in powers of H^, 
but the infinite diagmagnetic susceptiblity of a supercon- 
ductor replaces this with a non- analytic function of |. 
We will also discuss the case of k large but finite. 

Because the SC order is present on both sides of the 
transition, we can describe it in a static GL theory: the 
free energy Fq l is given by 



cPx 



iA] ijj 



X 



Here 'il'{x) is the complex SC order parameter, and we 
have performed standard rescalings to cast this the- 
ory in a dimensionless, universal form: all lengths (x) 
are measured in units of the (bare) superconducting 
coherence length ^, the vector potential is scaled as 
Wx X A = {H/H^2)^i energies are measured in units 
of i7^^^/(47r) [He is the field at which the free energy 
of the uniform superconductor equals that of the nor- 
mal state). When k ^ cxi, as is the situation in the 
experiments on the cuprates [Q, the screening of H by 
the supercurrents is negligible, and we can develop our 
theory for the SC to SC-I-SDW transition taking H to be 
uniform. We will not use any particular model for the dy- 
namics of Tp, apart from assuming that its slow evolution 
occurs on time scales longer than those of the spin fluc- 
tuations of interest, and is presumably associated with 
the thermal and quantum motion of vortices p^ . 



1 



"Normal" 




FIG. 1. Phase diagram of 5 + Fgl/T at T = and large 
N as function of the field H and the parameter r (which 
is similar, but not identical to the doping 5) in the limit 
K CO. This theory offers a complete, experimentally ap- 
plicable, description only of the SC to SC-I-SDW transition 
at small H, away from the tetra-critical point M. The re- 
maining phase diagram is qualitative, and the non-SC phases 
should have some additional charge ordering. The path of the 
experiments of |l] is denoted by the vertical arrow. The up- 
per critical field of the SC state clearly depends upon r, and 
i/°2 is its value at M. The point A is at the non-universal 
value r — Tc, but the remaining phase boundaries can be ex- 
pressed in terms of rc and the couphngs in S [^. The point 
M is at H/H°2 = I, r = + v, the boundary BM is at 
H/H^2 = 1 - -I- v{rc " r))/(4«), CM is at r = + u, and 
DM is at H/H°2 = 1 + AfwA/(87rc^), where A is the solution 
of -f (iV-u/(27rc^)A + V - r + rc = Q. Near M, the posi- 
tion of AM is given by H/Hc2 = 1 — "diirc + v — r) /v where 
i9i = 1.1596 - 0{v'^) = 1.1803 - 0{v'^)) for a triangular 
(square) Abrikosov fiux lattice (the 0{v'^) terms will be de- 
scribed elsewhere) . The small H behavior of AM is one of our 
main results, and is in nUM. 



The SDW order parameter is a = 3 component vec- 
tor, (jiaix, r), where a = 1 . . . N and r is imaginary time. 
The quantum fluctuations of the (j)a as a function of r 
are known to play an important role even in the insu- 
lating antiferromagnet, and must surely be included in 
the vicinity of a quantum transition in a low dimensional 
magnet. For simplicity, we wiU assume that (j)a is real, 
but our theory admits a simple generalization to incom- 
mensurate ordering transitions requiring a complex order 
parameter The dynamics of (f>a is described by the 
action [0 (in the same units as Fgl)'- 
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where the x, t dependence of 0q is implicit, '^{x) is r in- 
dependent as discussed above, and c, r, u, u are coupling 
constants, with < 4m. We have neglected fermionic 
excitations because momentum conservation constraints 
suppress their coupUng to 0q 11 | . The action is the same 



as that near SDW ordering transitions in insulators, with 
the simplest symmetry-allowed terms in powers of and 



spatial and temporal gradients. The tuning parameter 
r (which, presumably, increases monotonically with in- 
creasing doping, S) drives the theory from the SC-I-SDW 
phase (with {(f)^) ^ 0, (t/j) ^ 0) to the SC phase (with 
{(fia) = 0, (V') 7^ 0) with increasing r; at H = 0, T = 0, 
this transition is at r = Vc (Fig |l|) . The latter phase has 
a. S = 1 exciton (or a collective "spin resonance" ) associ- 
ated with oscillations of (j)a about (j)a — 0, with an energy 
which vanishes as r \ rc |jl^. Our main results for the 
H dependence of the physics in the vicinity of the SC to 
SC-I-SDW phase boundary depend crucially on the cou- 
pling V > between the amplitudes of the SC and SDW 
orders; such a repulsive coupling was emphasized in [^,^. 

A powerful tool to account for the quantum fluctua- 
tions of (pa is the large N expansion ||l^. At large N 
and small T, the saddle-point equations describing the 
properties of 5 -I- Fgl/T in phases with {(f>a) = are 
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-l + |0(a;)|2-(V, -iA)2j ij{x) 

+ {NvT/2) ^ G{x, X, ujn)i^{x) = 
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where ujn is a Matsubara frequency and G(a;, x' , uJn)Sap — 

fl/T 



dTe'^-^cj}aix,T)Mx',0)) is the 0^ 
tion which obeys 



Green's func- 



{luI- c^Vl + V{x))G{x,x' ,LUn) = S^ix^ x'). (4) 



The solution of (|2|j^,y) (and their straightforward gener- 
alization to phases with ^ |l^]) leads to the phase 
diagram in Fig |l|. We emphasize that present model is 
accurate only at small H in the vicinity of the SC to 
SC+SDW transition. The other phases in Fig |^ involve 
loss of SC order, and for these a more complete treat- 
ment of the charge fluctuations is surely needed: various 
site- and bond-centered charge ordered states ("stripes" 
and "spin-Peierls" ) and Wigner crystal states are likely 
to play a significant role (see e.g. |l^,|lj]) 

We now present an analytical description of the cru- 
cial physical properties of the solution of (||,§,^) in the 
SC phase at small H and close to a full numerical 
solution will be presented in future work. The SC or- 
der, ip, supports an Abrikosov flux lattice of vortices 
at TVu ^ GO positions {Ry} {il}{Rv) = 0), wit h a core 
size of order unity, and a lattice spacing Ly ~ V ^[^2/^- 
The resulting |-0(a;)p acts like a periodic potential for (pa, 
and it is useful express G in terms of the Bloch states, 
9nk{x) — Un{x)e'^'^^ of this periodic potential: 



G{x,x',u;n) = 

nk 



lki^)9nk{x') 



(5) 



where {-c^Vl +V{x))gnk{x) = elk9nk{x), n is a "band" 
index, k extends over the first Brillouin zone of the re- 
ciprocal lattice of {Rv}, and the eigenvalues e^^. > 0; 
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SDW order appears when one of e„fe first vanishes. An 
immediate experimental consequence of this structure is 
the presence of "Bragg reflections" at wavevectors sep- 
arated by the reciprocal lattice vectors in the dynamic 
spin structure factor, as shown in Figure 0. 
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FIG. 2. Dynamic spin susceptibility in the vortex state 

^nfc '^(^ ~ ^ frifc), integration is over the Bril- 

louin zone of the reciprocal lattice of {Kv\, n is a band index, 
Gm are reciprocal lattice vectors, and u„{Gm) are Fourier 
transforms of the periodic wavefunctions Note that 

dispersing S = 1 exciton modes S{(jj — enk) are shifted by all 
Gm with the weight \un{Gm)\^ ■ 

Our interest will primarily be in the nature of the low- 
est energy band, eofe: which also controls the transition to 
the SC-I-SDW state. As the cores of the vortices act as at- 
tractive potentials for gnk{x), a possibility for such a state 
is a superposition of Wannier orbitals localized at the 
cores of the vortices gok{x) = J2{r^.} e^^^" f{x — Rv) 
where f{x) is exponentially localized on a scale t <^ Ly. 
However, the self interactions of the (/>„, accounted for by 
the self-consistent potential proportional to u in (H), have 
a large impact on these localized states, and we will es- 
tablish that such a structure for gpfc must break down as 
the transition to the onset of SDW order is approached. 
The key argument was made by Bray and Moore [ pO[ in 
a different physical context, and we review their reason- 
ing. As we are assuming i <^ Lv, the localized states 
at the vortex cores can be treated independently of each 
other: eofc is independent of k and |/(0)| ^ Let 
V{x) = V> {V{x) — V<) for X outside (inside) the vor- 
tex cores where ip{x) ~ 1 {4'ix) ~ 0). Then subtracting 
evaluated at these x's from each other, and noting 
that the difference in the term proportional to u arises 
primarily from the localized state in the vortex core, we 
obtain (at T = 0) w - iVw|/(0)|Veoo ~ V> - V< > 0. 
This implies that the localization length must be at least 
as large as £ As eoo \ upon approaching 

the onset of SDW order, we must eventually have £ ~ L^, 
and the assumed structure for gok breaks down. 

We are, therefore, compelled to turn our attention to 
extended, plane-wave like wavefunctions for the gnk (we 
do not exclude non-trivial structure in these wavefunc- 



tions on the scale Ly). The self interaction terms are 
now expected to be less important: we initially neglect 
the terms proportional to G in (H,§) , but will account for 
them later. We now demonstrate that, for small H, the 
eigenenergies, eofc are controlled by universal structure in 
ipix) in the region 1 <C ^ well outside the vortex 
core at i?^ = (and in the corresponding regions of the 
other vortices). Analysis of (^), following the standard 
description of a superconducting vortex, shows that in 
such regions 

\^{x)\^l-l/{2x'); (6) 

the second term above is the correction in the amplitude 
of the SC order induced by the superflow kinetic energy. 
We now use perturbation theory to evaluate the change 
in eoo induced by (^); this requires the quantity 

mx)\'), = 1 {H/{2H^^)H{}^H^JH), (7) 

where the average is over spatial positions x, the log- 
arithm arises from the integral of l/x^ over the two- 
dimensional unit cell of the vortex lattice (it is cutoff 
at short scales by the superconducting coherence length 
~ 1, and at long scales by L^), and the constant z?2 ~ 3 
(for triangular and square vortex lattices) was obtained 
by numerical solution of (||) using the method of [^ . 
From (^), we obtain the leading perturbative correction 

^loiH) - 4(0) - {vH/{2H^^)) H^2H",,/Hy, (8) 

this result is also a variational upper bound on eoo(i?), 
associated with the wavefunction 500 — 1- For small 
r — Tc we have £00 (0) — ^ '"c) (this corresponds to the 
exponent v = 1/2 in mean field theory), and so we can 
rewrite (||) as eoo{H) — {r{H) — rc)^^'^ where 

r{H) ^ r - {vH/{2H^^)) H^^H^JH). (9) 

The vanishing of eoo{H) determines the position of the 
phase boundary AM between the SC and SC-I-SDW 
phases at small H: 

H/H",, = 2{r-r,)/[v\n{l/{r-n))]; (10) 

the variational argument shows that this is an upper 
bound for H, and so the phase boundary can only be- 
come flatter at higher orders. A fully self-consistent cal- 
culation which includes the terms proportional to G in 
(^,^ is more involved, but tractable: we find a modi- 
fied functional dependence of €oo{H) on r{H), eoo(-ff) = 
2TTc'^{r{H) — rc)/{Nu{l — v^/{4:u))) (corresponding to the 
exponent 1/ = 1 in the large N limit ||l8| ) , but the expres- 
sions (^) and ( p^ continue to hold. 

We comment further on experimental implications of 
our results. For very small if gofe is localized on a scale 
£ <ti Ly, the lowest energy magnetic mode is a relatively 
high energy S = 1 excitonic mode near the vortex cores. 
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as noted earlier. The spectral weight of this mode will be 
much smaller than the slightly higher energy bulk contri- 
bution from outside the vortex cores, and will consist of 
a feature of width Sk = (resolution of the Bragg re- 
flections may not be feasible). For slightly larger H, the 
distinction between the localized and delocalized contri- 
butions will disappear, and we expect that a full dynamic 
spin structure with "Bragg reflections" may be resolved 
in neutron scattering experiments, with the lowest energy 
mode obeying (||). It is also interesting to note that, even 
for small H, repulsion between the SC and SDW order 
parameters described by the v term in (|I|) will lead to an 
interesting structure in the spatial form of on the 

scale i which may be significantly larger than the bare 
SC coherence length (unity in our units) . This surprising 
effect, where the superconducting vortex core structure 
depends on the magnetic field for H ^ _ff°2, is another 
significant prediction of our theory, and is testable e.g. 
in tunnelling experiments. 

An interesting recently observed effect may read- 
ily be accounted for by 5 + Fcl/T: the intensity of 
the elastic neutron scattering increases when H is ap- 
plied to the SC-fSDW state (fine BA on Fig 0). In 
the large N limit we find that the expectation value 
of the staggered moment increases with H as (l^al)^ — 
(rc - r{H))/(2u{l - wV(4u))), with r{H) given in (|). 
Quantitative agreement with this relation has been re- 
cently observed Q. 

We now discuss the consequences of supercurrent 
screening of the magnetic field at finite k. For finite k, 
and small H we will have a Meissner phase where no vor- 
tices penetrate the sample. As shown in Fig ^, we may 
have Meissner phases that are purely superconducting or 
have coexisting magnetism. It is interesting to note that 
a finite density of vortices penetrates a sample at -ffci, so 
there is a finite jump in the SDW order across the lines 
BMi and M1M2. 




FIG. 3. Phase diagram of S + Fgl/T for finite k and small 
fields. The SC and SC-I-SDW phases in Fig |l| are identical to 
the corresponding "vortex " ph ases here. The position of Hci 
is given by Hci/H^2 = •\/7i"/2 (In K/K^)i/)|f^Q, where for the 
BMi line ipjj^o = 1 - v{rc - r)/{Au{l - v'^ /{Au))), and for the 
M1M2 line i/;|f^() = 1 -f Nv eoo(0)/87r with eoo discussed after 
equations (g) and (j^Bj). 

We conclude by noting some broader implications of 
our results, (i) In the discussion above we assumed that 
v>0, which is consistent with the experimental situation 



in La2-5Sr5Cu04. The case v <0, where the transition 
into the SC state leads to an enhancement in the SDW 
fluctuations, may be also described following the formal- 
ism above, (ii) Closely related phase diagrams should 
apply to other ordering transitions of superconductors 
in a field, including charge density wave and "staggered 
flux" orders. 
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